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Summary 

Computer-aided tomography (CAT) provides a 
means of noninvasively measuring the air-density 
distribution around an aerodynamic model. This 
technique is global in that a large portion of the 
flow field can be measured. A test of the applica- 
bility of CAT to transonic velocities was the 
objective of this investigation. A hemispherical- 
nose cylinder afterbody model was tested at a Mach 
number of 0.8 with a new laser holographic inter- 
ferometer at Ames Research Center's 2- by 2-Foot 
Transonic Wind Tunnel. Holograms of the flow field 
were taken and were reconstructed into interfero- 
grams. The fringe distribution (a measure of the 
local densities) was digitized for subsequent data 
reduction. A computer program based on the 
Fourier-transform technique was developed to con- 
vert the fringe distribution into three-dimensional 
densities around the model. Theoretical aerody- 
namic densities were calculated for evaluating and 
assessing the accuracy of the data obtained from 
the tomographic method. 

INTRODUCTION 

For many years, aerodynamicists have used 
Mach-Zehnder interferometry for both qualitative 
flow visualization and quantitative density mea- 
surements. However, that technique is not practi- 
cal in larger wind tunnels because of its sensitiv- 
ity to vibration and the expense of the very high 
quality optics it requires. Thus, its use has been 
limited to small laboratory experiments. 

During the past decade laser holographic 
interferometry has been developed into a practical 
technique that possesses excellent potential for 
large-field, quantitative density measurements. 

Many research laboratories are now using laser 
holographic interferometry for airfoil work, turbo- 
machinery research, jet-exhaust studies, and many 
other types of fluid-flow investigations [1-15]. 
These applications can be categorized into two- 
dimensional and three-dimensional flows. Two- 
dimensional flows would include airfoils [13,14], 
corner flows [15], and channel flows [ 1 6 ] . 

Holographic interferometry for two-dimensional 
flows has proved to be a practical measurement 
method, one reason being that the optical path 


lengths in two-dimensional flows are usually long 
enough so that there are sufficient fringes for 
good signal resolution. In addition, an algebraic 
expression that relates the density to the fringe 
count makes the data reduction problem manage- 
able. In fact, for the case of infinite fringe 
interferograms, the fringes are contours of con- 
stant density. 

For the general case of three-dimensional 
flow, such as airplane models, rotating rotor 
blades, Jets, and nozzles, CAT (Computer Aid Tomog- 
raphy) methods are being used to extract the den- 
sity from the fringe data. CAT combines holo- 
graphic interferometry, for taking multiple views 
of the flow, and a computer to solve an integral 
equation relating the fringe data to the density. 
Some initial experiments [12,17,18] indicate that 
the CAT methods are potentially effective means of 
measuring three-dimensional flows. Perry and Lee 
[19] tested a hemispherical-nose cylindrical after- 
body model at Mach number of 0.8 to assess a CAT 
method for axisymmetrical flow. The emphasis in 
that work was on the experimental method, recon- 
struction of the holograms, identification and 
digitization of the fringe pattern, and a Fourier- 
transform algorithm to obtain density from the 
fringe data. 

This paper presents the analysis of the com- 
plete flow field and compares the results with 
theoretical calculations to assess the usefulness 
of this CAT method. 

EXPERIMENTAL METHOD 

Interferometer 

The holographic interferometer is an off-axis, 
two-beam system. It was specially constructed in 
1981 for permanent installation in the 2- by 2-Foot 
Transonic Wind Tunnel at Ames Research Center. The 
interferometer uses a pulsed Nd:yag laser to pro- 
vide a 532-nm output. It was modified by enclosing 
the object beam in a cloth bellows for the purpose 
of reducing room air currents, which produced noisy 
fringes, thereby yielding poorer data. The refer- 
ence beam was enclosed in a concrete trench. A 
complete description of the interferometer is given 
by Craig et al . [ 10] . 
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Model 

The model was a stainless steel hemisphere 
cylinder, 5.08 cm in diameter and 35.56 cm long. 

The model was mounted on the tunnel sting support 
and set to within 0.01° of the centerline of the 
test section in order to obtain axisymmetric 
flow. With axisymmetric flow, the density field 
about the model axis is symmetrical and only one 
hologram taken normal to the flow is necessary to 
define the flow. 

Holographic Method 

Dual-plate interferometry [20] was used to 
record the holograms. Several reference holograms 
were taken at zero velocity and several object 
holograms were taken at the test velocity. From 
these, at least one set having sufficient fringe 
contrast, brightness, and sharpness was obtained. 
The holograms were reconstructed into infinite- 
fringe interferograms (Fig. 1). In interpreting 
the interferograms, zero was assigned to the large, 
bright fringe either in the free stream or in the 
area well ahead of the model. All subsequent 
bright fringes were assigned the numbers 1, 2, 

3, ..., consecutively. The semiautomatic method of 
Becker and Meier [21], which relies on an image 
processor and a minicomputer, was used to digitize 
the fringes on planes normal to the model axis 
(Fig. 2a). Van Houten's [22] Four ier-transform 
method was used to obtain the density by inverting 
an integral equation relating the index of refrac- 
tion (hence the density) to the fringe shifts. The 
method requires equally spaced views of the model 
over an arc of 180°. According to Sweeney [23], 
the accuracy increases with more views but so does 
the computer time. It was decided that 72 views 
would be a good compromise. (Remember that for 
axisymmetric flow, only one view is needed to 
define the flow. This view was repeated 72 times 
in the Fourier method.) 

RESULTS 

Figure 1 shows a typical interferogram used 
for data reduction. Examination of the fringe 
pattern shows some asymmetry in the flow; it was 
typical of many obtained at this test condition. 

The major density changes occurred within a radius 
of 2 model diameters from the model nose tip. 
Therefore, the fringes were digitized on 13 planes 
in this region. The planes were spaced equally 
one-quarter model diameter apart and were normal to 
the model axis, as shown in Fig. 2b. Since the 
model was set for axisymmetric flow, the fringe 
data were averaged over the upper and lower half of 
the interferogram for inversion purposes. The 
apparent asymmetry was due to noise caused by a 
combination of laser instabilities, reconstruction 
sensitivity, wind-tunnel vibration, flow 
angularities, beam wandering, and unsteady shock 
waves on the model . 


The most important noise source was the laser 
beam wandering. This wandering occurs over a period 
of several minutes or longer. For the worst case, 
the beam can move so much that it will miss half of 
the hologram and make the hologram unusable. Since 
it takes the tunnel several minutes to reach test 
condition, the degree of beam wandering within this 
time period will determine the degree of noise. 

Some noise can be seen in the region ahead of the 
model and is exhibited by broad random fringes 
(Fig. 1). The order of the noise is one fringe and 
it exists even for ambient flow. This means that 
the zero fringe, required for accurate results, 
cannot be determined to within a resolution of one 
fringe. To resolve this problem, this first inver- 
sion plane has to cross a point where the density 
is known. In this experiment, the known density 
was at the model nose-tip. This density is the 
isentropic density at the stagnation condition. 

The fringe distribution at the nose-tip plane 
is shown in Fig. 3. By inverting this fringe dis- 
tribution and matching the density at the nose tip 
with the known density, the zero fringe was 
found. This was done by an iterative technique and 
was accomplished in two or three iterations. Once 
the zero fringe was found from the stagnation-point 
solution, fringes on the entire interferogram were 
identified. (It should be noted that for a model 
for which the stagnation point is not easily iden- 
tified, a known density at some point on the model 
could be used to determine the zero fringe.) 

The inversion process in the region ahead of 
the nose was straightforward. Behind the nose, 
fictitious fringes were added in the region where 
the model blocks the beam because the inversion 
theory requires continuous data over the entire 
plane. The inversion code was written with a 
spline curve-fitting routine to fill in the ficti- 
tious fringe data. According to Zien et al. [8], 
the density solution was not affected by the ficti- 
tious fringe distribution for axisymmetric flow. 

This assumption was checked by assuming three radi- 
cally different fictitious fringe distributions. 

The fictitious data from the spline method are 
plotted in Fig. 4a. The other two curves were 
chosen to be radically different from the spline 
curve. The resulting density distributions are 
shown in Fig. 4b. The densities in the region not 
blocked by the model were not affected by the dis- 
tribution of fictitious fringes assumed in the 
shadow of the model. 

Comparisons of the density data with theoreti- 
cal results are shown in Figs. 5a through 5i. Two 
computations, those of Coakley [24] and Pulliam 
[25], were used to assess the accuracy of the 
data. Both methods solved the steady, inviscid 
Euler equations. Pulliam's method can handle 
three-dimensional flows and was used in anticipa- 
tion of non-ax isymmetric cases. Coakley 's method 
can handle only axisymmetric flows; however, it 
needs only to solve the flow on one plane so, for 
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equivalent computer time, the mesh size can be 
smaller and the solution more accurate. Coakley's 
mesh size was set to be nearly twice as dense as 
that of Pulliam's. 

The comparisons of the experimental and theo- 
retical results in the region ahead of the model 
are presented in Fig. 5a to 5d; as shown in those 
figures, the experimental data are within \% of the 
theoretical results. Note that there was some 
scatter in Pulliam's results, a result of the 
coarse mesh he used. Referring back to Fig. 1, the 
signal generated in the region ahead of the model 
was about four fringes, as characterized by the 
four nearly elliptically shaped fringe patterns. 
Becker's automatic method digitized the line seg- 
ment on each bright fringe and used a spline fit to 
form a smooth curve through the data points. A 
maximum of 512 pixels or points to curve fit are 
available in the system. As the results show, very 
good accuracy was obtained at a signal of four 
fringes and a noise of one about one fringe. How- 
ever, at lower signal levels, the accuracy should 
deteriorate. Lower signals will occur at reduced 
Mach and Reynolds numbers (based on lower tunnel 
pressures), and in regions of low density 
gradients. 

An attempt was made to provide more fringe 
data by reconstructing the interferogram into the 
finite-fringe mode, as shown in Fig. 6. By super- 
imposing a set of parallel "wedge" fringes on the 
fringe pattern, more fringe data points were 
obtained. This method in effect provided an inter- 
polation between the fringes. More parallel 
fringes can be obtained by increasing the angle 
between the holograms during reconstruction. At 
some point, adjacent fringes will tend to merge 
because of the noise. The interferogram shown in 
Fig. 6 can probably hold twice as many fringes 
before merging occurs. The noise can be seen by 
the degree of waviness of a fringe from a straight 
line. Without noise, the fringes in the free 
stream would be straight lines. The noise was 
found to vary by as much as a fringe. Digitizing 
the finite-fringe data was a tedious process 
because the automated scheme was not designed for 
that purpose. The results of this method did not 
seem to improve the fringe data over those of the 
other method and it was subsequently dropped. 

Comparison of the experimental data and theo- 
retical results at the nose-tip plane is given in 
Fig. 5e; the agreement was excellent (within \l ) . 

At the nose tip, or stagnation point, Pulliam's 
result was 2% low because of the coarse mesh he 
used. Note the agreement with Coakley's result and 
remember that the solution was normalized to the 
isentropic stagnation-point density. 

Figures 5 present the comparisons on the hemi- 
sphere-nose (Z/R = 0.5, Fig. 5f; and 1.0, Fig. 5g) 
and just downstream of the hemisphere-nose 
(Z/R s 1.5, Fig. 5h). The agreement between 
experiment and theory was poor in a region within 


2 nose radii of the model. This poor agreement was 
due to a system of unsteady weak shocks and local 
flow separation which could not be predicted by a 
steady, inviscid theory. The shocks, which 
appeared as wrinkles on the fringes, were very 
difficult to see and were almost masked by noise, 
which also caused wrinkled fringes. 

In order to see the shocks better, a shadow- 
graph was reconstructed from the hologram 
(Fig. 7). The shocks and the flow separation can 
now be readily seen. A lambda shock was centered 
at about Z/R = 1.3 and the shock on the upper 
half of the model seems stronger. The boundary- 
layer separation started at about Z/R - 0.8. 

These results were also observed by Hsieh [26] at 
Mach numbers of 0.7 to 0.9. Coakley and Pulliam 
predicted single shock waves located at distances 
of 1.3 and 1.2, respectively. The flow was 
observed during the experiment by operating the 
interferometer as a shadowgraph at 10 to 20 pulses 
per second. The flow was seen to be unsteady and 
the shocks oscillated and varied in strength with 
time. The data in Figs. 5f-5g showed that beyond 
2 nose radii, or beyond the influence of the 
shocks, the agreement between experiment and theory 
was again within 

Figure 5i shows a typical density profile in 
the cylindrical portion of the model. The flow in 
this region approaches uniform free-stream den- 
sity. The data approached the free-stream values 
faster than the predictions. This is probably a 
result of the deterioration of the experimental 
accuracy. The signal-to-noise ratio in the region 
is about 1. The maximum deviation between the data 
and the predictions was about 2%. 

CONCLUDING REMARKS 

This computer-aided tomography experiment, 
conducted in a large wind tunnel using a state-of- 
the-art laser holographic system and a recently 
developed semiautomated digitizing system, demon- 
strated the feasibility of using tomographic 
methods to make accurate density measurements. 

Three changes are indicated that would improve the 
quality of the holograms, make the reconstruction 
easier, and speed up the data process. First, a 
laser with improved stability would reduce the beam 
wandering and noise and improve accuracy. Second, 
the use of newly developing thermoplastic recording 
devices would permit faster and more accurate 
reconstructions. And third, a positive method for 
reconstructing infinite-fringe interferograms would 
also improve accuracy. 

ACKNOWLEDGMENTS 

The author would like to thank Dr. R. L. 

Perry, University of the Pacific, Stockton, Calif., 
for his assistance in reconstructing the holo- 
grams. Drs. T. J. Coakley and T. H. Pulliam of 


3 



Ames Research Center were very helpful in setting 
up their numerical codes and gave valuable consul- 
tations. Dr. F. Becker, NRC, Ames Research Center, 
wrote the automated digitizing routines for this 
experiment and assisted in the initial 
digitization . 

REFERENCES 

1. Bemish, J. K. et al., "Wind-Tunnel Diag- 
nostics by Holographic Interferometry," AIAA Jour- 
nal , Vol. 7, Oct. 1969, pp. 2041-2042. 

2. Matulke, R. D. and Collins, D. J., "Deter- 
mination of Three-Dimensional Density Fields from 
Holographic Interferometry," Journal of Applied 
Physics , Vol. 42, Mar. 1971, pp. 1109-1119. 

3. Jogota, R. C. and Collins, D. L., "Finite 
Fringe Holographic Interferometry Applied to a 
Right Circular Cone at Angle of Attack," Journal of 
Applied Mechanics , Vol. 39, Dec. 1972, pp. 897-903. 

4. Reinheimer, C. L., "Holographic Subsonic 
Flow Visualization," Applied Optics , Vol. 9, Sept. 
1970, pp. 2059-2065. 

5. Hirth, A., Smigielski, P., and Stimpfling, 
A., "Use of Holography for Visualization of the 
Wake of Projectiles in Hypersonic Flight at 

Mach 6," Optics and Laser Technology , Nov. 1971, 
pp. 195-199. 

6. Redley, R. J., Jr. and Havener, A. G., 
"Application of Dual Hologram Interferometry to 
Wind-Tunnel Testing," AIAA Journal , Vol. 11, Nov. 
1973, pp. 1332-1333. 

7. Belozerov, A. F. et al., "Hologram Study 
of Gas Flow in a Ballistic Wind Tunnel," Soviet 
Physics - Technical Physics , Vol. 18, Oct. 1973, 
pp. 488-489. 

8. Zien, T. F., Ragsdale, W. C., and Spring, 

W. C., "Quantitative Determination of Three- 
Dimensional Density Field by Holographic Interfer- 
ometry," AIAA Journal , Vol. 13, July 1975, 

pp. 841-842. 

9. Surget, J., "Holographic Interferometer 
for Aerodynamic Flow Analysis," International Sym- 
posium of Flow Visualization, Bochum (RFA), Hemi- 
sphere Pub. Co., Sept. 1980. 

10. Craig, J. E., Lee, G., and Bachalo, 

W. D., "Nd:yag Holographic Interferometer for 
Aerodynamic Research," SPIE Paper 353-15, San 
Diego, Calif. , Aug. 1982. 

11. Burner, A. W. and Goad, W. K., "Flow 
Visualization in a Cryogenic Wind Tunnel Using 
Holography," NASA TM-84556, 1982. 

12. Matulka, R. D., "The Application of Holo- 
graphic Interferometry to the Determination of 
Asymmetric Three-Dimensional Density Fields in Free 
Jet Flow," Ph.D. Thesis, U.S. Naval Postgraduate 
School, Monterey, Calif., June 1970. 

13. Johnson, D. A. and Bachalo, W. D., "Tran- 
sonic Flow Past a Symmetric Airfoil--Inviscid and 
Turbulent Flow Properties," AIAA Journal , Vol. 18, 
Jan. 1980, pp. 16-24. 


14. Spaid, F. W. and Bachalo, W. D., "Experi- 
ments on the Flow About a Supercritical Airfoil, 
Including Holographic Interferometry," AIAA Paper 
80-0343, Pasadena, Calif., 1980. 

15. Kosakoski, -R. A. and Collins, D, L., 
"Application of Holographic Interferometry to Den- 
sity Field Determination in Transonic Corner Flow," 
AIAA Journal , Vol. 12, June 1974, pp. 767-700. 

16. Delery, J. M. , "Investigation of Strong 
Shock Turbulent Boundary Layer Interaction in 2-D 
Transonic Flows with Emphasis on Turbulence Phenom- 
ena," AIAA Paper 81-1245, Palo Alto, Calif., 1981. 

17. Clark, L. T. et al . , "A Three-Dimensional 
Density Field Measurement of Transonic Flow from a 
Square Nozzle Using Holographic Interferometry," 
Journal of Fluid Engineering , Vol. 99, Dec. 1977, 
pp. 737-744. 

18. Kittleson, J. K., "A Holographic Inter- 
ferometry Technique for Measuring Transonic Flow 
near a Rotor Blade," 9th European Rotorcraft Forum, 
Stresa, Italy, Sept. 13-15, 1983. 

19. Perry, R. L. and Lee, G., "Interferometry 
Applied to Symmetric Aerodynamic Models in a Wind 
Tunnel," SPIE Conference on Holography, 

Los Angeles, Calif., Jan. 1985. 

20. Vest, C., Holographic Interferometry , J. 
Wiley & Sons, New York, 1979- 

21. Becker, F. and Meier, G. E. A., "Entwick- 
lung eins Gerates zur Auswertung von Interferogram- 
men," Max-Planck-Institut fur Stromungsforsehung , 
Bericht 1 10/1981, June 1981. 

22. Van Houten, P. E., "The Application of 
Holographic Interferometry to the Determination of 
Discontinuous Three-Dimensional Density Fields," 
Thesis for Aeronautical Engineer Degree, Naval 
Postgraduate School, Monterey, Calif., 1972. 

23. Sweeney, D. W. , "Interferometric Measure- 
ment of Three-Dimensional Temperature Fields," 

Ph.D. Thesis, U. of Michigan, Ann Arbor, Mich., 
1972. 

24. Coakley, T. J., "Implicit Upwind Method 
for the Compressible Navier-Stokes Equations," AIAA 
Journal . Vol. 23, No. 3, pp. 374-380, 1985. 

25. Pulliam, T. H., "Implicit Finite Differ- 
ence Simulations of Three-Dimensional Compressible 
Flows," AIAA Journal , Vol. 18, No. 2, pp. 159-167, 
1980. 

26. Hsieh, T., "An Investigation of Separated 
Flows about a Hemisphere-Cylinder at Incidence in 
the Mach Number Range from 0.6 to 1.5," AIAA Paper 
77-179, Los Angeles, Calif., 1977. 

GEORGE LEE received his bachelor's and master's 
degrees in mechanical engineering from the Univer- 
sity of California at Berkeley in 1954. He has 
been at NASA Ames Research Center since that 
time. His principal research interests have been 
in aerodynamics, lasers, and interferometry. 


4 





1.27 cm 


PLANES EQUALLY SPACED 
1/2 NOSE RADIUS 


b) 13 planes used for data reduction covering 

Z/R = ±3. 


Fig. 2 Planes used for digitization. 
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Fig. 4 Effect of fictitious fringes in the region 
of the model, a) Fictitious fringe distributions; 
b) density ratio. 
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Fig. 6 A typical finite-fringe interferogram. 


Fig. 7 Shadowgraph showing shock waves. 
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